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Persistent spin dynamics intrinsic to amplitude-modulated long-range magnetic order 
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An incommensurate elliptical helical magnetic structure in the frustrated coupled-spin-chain sys- 
tem FeTe20sBr is surprisingly found to persist down to 53(3) mK (T/Tat ~ 1/200), according to 
neutron scattering and muon spin relaxation. In this state, finite spin fluctuations at T — > are 
evidenced by muon depolarization, which is in agreement with specific-heat data indicating the 
presence of both gapless and gapped excitations. We thus show that the amplitude-modulated mag- 
netic order intrinsically accommodates contradictory persistent spin dynamics and long-range order 
and can serve as a model structure to investigate their coexistence. 

PACS numbers: 75.25.-j, 75.30.-m, 75.30.Fv, 75.40.Gb, 76.75.+i 



u 

-f— > 

a 

1 

C 

O 

o 



> 

00 
in 

00 
o 



Geometrical frustration is a common precursor for ex- 
otic magnetic ground states - from long-range ordered 
(LRO) incommensurate (IC) spiral states to highly dis- 
ordered frozen spin-glass states It can even lead to 
spin liquids, where spin fluctuations may endure down to 
zero temperature - persistent spin dynamics (PSD) Q 
is present. Such behavior is typically found in highly- 
frustrated pyrochlore and kagome spin systems with 
macroscopically degenerate ground states [3j . Since spin 
fluctuations hinder the onset of extended static correla- 
tions, PSD and LRO are generally considered mutually 
exclusive. Remarkably, their coexistence in the same 
phase has been reported in several frustrated magnetic 
systems [H-fl], but still lacks a suitable explanation. 

To explore this phenomenon we focus on the 
FeTe205Br multiferroic, in which the magnetic exchange 
network consists of alternating Fc 3+ (S — 5/2) spin chains 
coupled by weaker frustrated interactions within the be 
layers @ (Fig.[l|). The magnetic order at 5K, i.e., well 
below the Neel temperature Tjy = 10.5 K, was described 
as an IC collinear amplitude-modulated (AMOD) struc- 
ture, with magnetic vector q=(| 0.4630) In such 
state only part of the total Fe magnetic moment at each 
site contributes to LRO, while its counterpart is disor- 
dered and may fluctuate. Generally, on cooling the or- 
dered component fully develops, which manifests either 
as a " squaring" of the AMOD structure llj or as a con- 
secutive transition, where a perpendicular ordered com- 
ponent develops leading to a circular helix [l2j . Since 
FcT^OsBr exhibits no subsequent transition down to 
300 mK 13(, it is possible that PSD and LRO coexist. 



In this letter, we report on a combined study of spher- 
ical neutron polarimetry (SNP) and neutron diffrac- 
tion, which reveals that the IC AMOD magnetic struc- 
ture persists to the lowest accessible temperatures 
(T/T/v ~ 1/175). This is consistent with muon spin re- 



laxation (/iSR) measurements at T/Tn ~ 1/200, which in 
addition to static LRO signify the presence of PSD. The 
coexistence of LRO and PSD is supported by specific- 
heat data, indicating gaped as well as gapless magnetic 
excitations. Our study suggests that this is intrinsic 
to AMOD magnetic structures. It offers a well-defined 
framework and a coherent explanation for the coexistence 
of LRO and PSD that has been missing for the known 
cases [8|. 
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FIG. 1: Top: Magnetic structure in FeTe 2 5 Br at 2K. Or- 
ange and blue spheres are Fei and Fe 2 magnetic ions, re- 
spectively. The color of the magnetic moments (arrows) are 
different for the eight magnetically inequivalent sites Fe m „. 
The amplitude of the moment changes with its orientation 
within the ellipsoidal envelope. Black and gray lines denote 
the dominant intrachain and the weaker interchain exchange 
interactions, respectively. Bottom: IC amplitude modulation 
of the Fen magnetic moment along the b-axis. 



TABLE I: Parameters of the best magnetic structure model at 
2K for two independent magnetic atoms (Fei and Fe2), and 
eight magnetic phases ipmn, i-e., one for each of the magnetic 
Fe mn atoms in the unit cell (m = 1,2, n = 1-4). The sites Fei 2 - 
Few are obtained from Fen [0.1184(6), -0.001(1), 0.9734(7)] 
and Fe 2 2-Fe 2 4 from Fe 2 i [0.9377(6), 0.2953(1), 0.8562(6)] by 
symmetry elements i, 2i„ and 2i y i, respectively. The orien- 
tation of the moments is given in the a* be coordinate system. 
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We begin by a re-investigation of the FeTe205Br 
magnetic ground state, combining single-crystal neutron 
diffraction and SNP, which was proven before to be in- 
valuable for determination of complex magnetic struc- 
tures The experiments were performed at 5 and 
2K, respectively, on high-quality single crystals [lj with 
an average size of 15 x 8 x 2 mm 3 at the Swiss Neu- 
tron Spallation Source (SINQ), Paul Scherrer Institute 
(PSI), Switzerland. For SNP the MuPAD device on 
the triple axis spectrometer TASP (A = 3.2 A) was used, 
while the diffraction experiment employed the single- 
crystal diffractometer TriCS (A = 2. 317 A). Polarization 
matrices were measured for three different crystal orien- 
tations [IH, accessing for the first time [3, EH also hkl ± q 
(1 ^ 0) magnetic reflections. 

The general magnetic structure model, in which the 
magnetic order breaks all crystallographic symmetry op- 
erations [l|, [l6[ , dictates the magnetic moment at a par- 
ticular Fe site to follow an elliptical helix with the pitch 
along the magnetic q vector, 



S m „(r;) = Somn cos(q-r 4 — ij} n 



)+Sq m 



l(q-ri-Vw)- (1) 



Here, vector r.; defines the origin of the z-th cell, m = 1,2 
identifies the crystallographically inequivalent Fe-sites, 
and n=l-4 denotes the four Fe positions within the crys- 
tallographic unit cell (see caption of TableU]). The com- 
plex vector So mn is determined by its real and imag- 
inary components, and So^ n , defining the am- 
plitude and the orientation of the magnetic moments, 
while ipmn denotes a phase shift. We assume the same 
moment So mn = So m for all crystallographically equiv- 
alent Fe-sites. The reliability of the refined magnetic 
structure is ensured by simultaneous refinement of all 
SNP and integrated magnetic-peak intensities data flEl ]. 
The best solution (Fig. [T] and TableHJ with the dominant 
component pointing ~ 45 ° away from a towards -b axis 
agrees with the previously proposed simplified collinear 
IC AMOD structure In addition, a small perpen- 
dicular component is identified, resulting in an overall 
elongated-elliptical cycloid (|S^°|/|S m | - 0.37), with its 
normal canted ~ 15 ° away from the c-axis. The new data 
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FIG. 2: Refinement quality for (a) polarization matrices and 
(b) integrated intensities, (c) Top: fc-scans at 60(3) mK and 
2K across (3.5 -0.5±<5 0) reflections, witnessing the absence 
of higher harmonics, expected to occur at marked positions. 
Bottom: the difference between the two scans which match 
within the error bar. 



thus reveal a magnetic structure that combines AMOD 
and helical properties. The quality of the new refine- 
ment reflects in the reduced x 2 = 9.6 of the polariza- 
tion matrices being notably decreased with respect to its 
value in the collinear AMOD (14.7) and the circular he- 
lical structures (14.9) A small misfit [Figs.[5Ja) and 
(b)] inevitably originates from a weak nuclear-magnetic- 
interference term |15j . 

At T — > 0, the ordered component of the magnetic mo- 
ments in the elliptical IC AMOD structure is expected 
to grow at the expense of the disordered one. This can 
reflect either as a change of the magnetic-reflection inten- 
sities, if circular helix is formed, or alternatively, as addi- 
tional magnetic reflections with propagation vectors 3q, 
5q, ... and intensities 1/9, 1/25, ... of the first-order mag- 
netic reflections in case of " squaring" [ijj . To detect such 
changes, single crystal neutron diffraction was performed 
at 60(3) mK {T/T N ~ 1/175), where we recorded broad 
fc-scans of several hkl ±q pairs, i.e., intersecting the po- 
sitions of hkl ± nq, n = 3, 5, ... reflections. Surprisingly, 
we find no significant difference between 60 mK and 2 K 
data as well as no trace of higher harmonics [Fig.^c)] , 
which implies that both the AMOD ordered component 
and its disordered counterpart are still present at the 
lowest accessible temperatures. 

Since the neutron scattering experiments could probe 
only static magnetism, we employed a local-probe /xSR 
technique, which is extremely sensitive to internal mag- 
netic fields and can distinguish between fluctuating and 
static magnetism, as well as between LRO and static 
magnetic disorder Q. The /iSR-experiments were per- 
formed on the MUSR instrument at the ISIS facility, 
Rutherford Applcton Laboratory, United Kingdom, on 
the same high-purity powder samples as used in our ear- 
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lier study [15J, ll8| . All the data in the following is shown 
with properly subtracted background signal (~ 15 %) 
fisT ]. Preliminary results indicated that muons stop at 
several inequivalcnt positions [18| . However, to test the 
magnetic structure model, these must be precisely de- 
termined. New measurements were therefore first per- 
formed in the paramagnetic state, at 50 K [Fig.[3ta)], 
where weak static nuclear magnetic fields are expected 
to govern the fi + spin relaxation Q. Since these fields 
can be exactly calculated from the crystal structure, 
they are essential for identification of the muon stop- 
ping sites, as demonstrated below. In case of a single 
muon stopping site, in a paramagnetic powder sample 
the muon polarization in zero applied magnetic field (ZF) 
is given by the Gaussian Kubo-Toyabe relaxation func- 
tion G KT (i,A) = i + |[1 - (Ai) 2 ]oxp[-(Ai) 2 /2] multi- 
plied by the exponential function exp[— Aj,i]. The former 
accounts for a static Gaussian nuclear field distribution 
with the width AAy^ (7^ = 2ttx 135.5 MHz/T), whereas 
the latter describes a weak dynamical relaxation due to 
fast electronic fluctuations [8[- The resulting function 
has a single dip, which is removed only by a fast dy- 
namical relaxation, i.e., Al>A. Our ZF data exhibits 
a more complex behavior [Fig.(3{a)], as they require two 
two-component model 

G{t) = [G KT (t, A A ) + G KT (i, A B )] • exp[-Ait]. (2) 

This model is supported by measurements in longitu- 
dinal applied magnetic fields (LF), where decoupling 
of the muon relaxation from nuclear-magnetic fields 
occurs in two steps at ^0.5 and ~2mT [Fig.^a)], 
respectively. Simultaneous fit of all 50 K /iSR data 
to the two-component model, extended for the case 
of LF [f| [Fig.|3fa)], yields the relative occupan- 
cies of the two muon stopping sites of 83(3) % and 
17(3)% with corresponding Aa/j^ = 0.064(5) mT and 
A B /7 M = 0.304(5) mT, whereas \ L = 0.039(3) /xs" 1 . 

Since /i + are positively charged particles and are gen- 
erally expected to stop at the electrostatic-potential min- 
ima, we calculated the electrostatic potential in the 
FeTe205Br unit cell using density functional theory fljj ]. 
This way several possible stopping sites were identified 
[Fig. [2th)], for which dipolar nuclear-magnetic field dis- 
tributions were calculated [l5j]. The dominant Aa/^^ is 
found to be in excellent agreement with nuclear field dis- 
tributions calculated at two local electrostatic-potential 
minima, Pi = (0.59, 0.12, 0.43) and P 2 = (0.25, 0.15, 0.90) 
[Fig-E^b)]. These are thus assigned as the prime muon 
stopping sites. On the contrary, Ae/j^, docs not 
agree with calculated distribution at any electrostatic- 
potential-minima [15j , so the third (least occupied) muon 
stopping site remains unassigned. 

Identification of the Pi and P2 stopping sites allows 
us to calculate local dipolar magnetic fields from the or- 
dered Fe moments and thus to double-check the magnetic 
structure determined by neutrons. We computed (l5| 




FIG. 3: (a) Longitudinal-field /^SR measurements on pow- 
der sample FeTe2 0sBr at 50 K (symbols) and corresponding 
fits (lines) - see text, (b) Electrostatic potential at 6 = 0.15, 
i.e., intersecting the most pronounced minima, including Pi 
and P2. (c) Zero-field fiSH measurement at 1.6 K with corre- 
sponding two- (dashed line) and three-component (solid line) 
fits given by the magnetic structure (see text). Insets: cal- 
culated field distributions Dij(B) at Pi (i = 1,2) sites for the 
four magnetically inequivalent positions j within the unit cell. 



normalized field distributions Dij(B) at four inequiva- 
lent sites (j = 1-4) [Wyckoff position 4(e)] for each P, 
(£=1,2) felt by muons stopping in random unit cells [in- 

^ is 

line) to the corresponding model 



sets in Fig.[3fc)1. In Fig.^c), the ZF /J.SR data collected 
at 1.6 K at PSI LL8fl is shown together with the fit (dashed 



G(t) = ±cxp[-(A L i) Q ] + |^[^exp(-Aii) 
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(3) 



In this expression the first term, commonly called the 
"1/ 3- tail" |8( , describes muons in a powder sample, whose 
initial polarization is parallel to the internal magnetic 
field and therefore changes only due to fluctuations of 
this field. This term (dynamical relaxation rate A/, and 
stretch exponent a) is determined from long-time decay 
measurements presented in the next paragraph. Now 
we focus on the second term, which depicts oscillations 
due to internal fields induced by the LRO magnetic or- 
der. These oscillations are damped (Ai) by spin fluctua- 
tions and/or static relaxation resulting from a distribu- 
tion of muon stopping positions around Pi [IB]], with oc- 
cupancy A4 Ai = 1). The fit yields Ai = 0.60(5) /txs" 1 , 
A 2 = 50(8) ^s" 1 and A2/A1 = 1.4(1). Most importantly, 
our model excellently accounts for the experimental os- 
cillation frequencies of the muon polarization determined 
by Dij(B) (no adjustable parameters) and thus firmly af- 
firms the elongated-elliptical IC AMOD state. The small 
discrepancy [Fig.[3{c)] is most likely due to the neglected 
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17(3) % of muons, with unknown stopping position (As). 
Indeed, adding a third component with A3 — 17% and 
its decay approximated by exp(— X 3 t) leads to a per- 
fect agreement with the experiment [Fig.[3Jc)]. The im- 
proved model yields A 3 = 11(1) /is -1 , while parameters 
Ai, A 2 and A 2 /Ai stay within the error bars of the two- 
component model. The agreement of the neutron and the 
/iSR results proves that ^tSR probes the intrinsic mag- 
netic properties and can thus provide an invaluable in- 
sight into spin dynamics of the FeTe2C>5Br system. 

In case of completely static (X L =Q) local magnetic 
fields, the "1/3 tail" persists in the t— >oo limit [Eg. (EH. 
Since preliminary ZF /iSR data implied its decay 18| , we 
extended these measurements to longer times and very 
low temperatures [Fig.@|a)]. The new data confirm the 
decay of the "1/3 tail" and clearly show its persistence 
down to the lowest accessible temperature of 53(3) mK, 
i.e., T/Tjv ~ 1/200. Since this decay can only be of dy- 
namical origin, it unambiguously proves that muons, ex- 
periencing static magnetic fields due to the IC AMOD 
order, experience also local- field fluctuations, i.e., reveal- 
ing the coexistence of PSD and LRO at T — > 0. 

To obtain deeper insight into PSD, we focus 
on the dynamic part of the /iSR signal, i.e., 
muon polarization for t > 0.5 /xs, where the sec- 
ond term in Eq. ([3]) is already relaxed [Fig.EJb)]. 
These data are fitted with G(t) = \ cxp[— (\ L t) a ] 
[Fig. Eta)], where a accounts for the dynamical re- 
laxation rate distribution p\ L (v), related to the 
stretchcd-exponcntial function by the La plac e transform 
exp[-(A L t)°] = p XL (is)atp[-(vt)}dv A small 

a = 0.30(2), found temperature independent below TV, 
indicates a very broad p\ L {v) [inset in Fig. Eta)], which 
most likely reflects a distribution and dynamical nature 
of the disordered parts of the magnetic moments in the 
AMOD state. Additionally, p\ L (v) can be broadened 
because three different muon stopping sites are present. 
The obtained Xl [Fig.HJb)] shows a linear temperature 
dependence and, most importantly, converges to a finite 
zero-temperature value A° = 0.012(2) /is -1 , characteristic 
of PSD. We note that in spite of the small A°, the broad 
p\ L (y) spans far into the pSR time window (v >0.1 /is -1 ), 
as evident from substantial experimental decay of the 
"1/3-tail" at early times [Fig.gfa)]. 

In LRO antifcrromagncts, \l(T) is expected to fol- 
low either T n dependence with n > 2 for T^>e g or 
cxp(-e ff /T) dependence for T<e ff [1 [U-Elj] . Here e g is 
the magnon energy gap. Considering e g — 11.5 K, as de- 
termined in the recent antifcrromagnctic-resonance study 
Q , Al (T) should change exponentially with temperature 
in the inspected temperature range, which is clearly not 
the case [Fig.[4|b)]. Hence, Al(T) signifies magnetic exci- 
tations, which are different from the usual magnon modes 
in LRO states. Finite relaxation at T^O is typically 
found in gaplcss spin-liquids where it is ascribed to quan- 
tum spin fluctuations @, H, El, - In FeTe 2 5 Br, these 
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FIG. 4: (a) Muon depolarization (symbols) at low tempera- 
tures, indicating the presence of PSD, and corresponding fits 
(lines) to the stretched-exponential function exp[— (Al£) ]• 
Inset: relaxation rate distribution at 53 mK determined from 
fit (see text), (b) Extracted relaxation rates, showing linear 
temperature dependence (solid line). For comparison activa- 
tion behavior (dashed line) for gapped (e 9 = 11.5 K) magnon 
excitations is shown, (c) Temperature dependence of the mag- 
netic contribution to the specific heat, c mag (T), together with 
corresponding fits (see text for details). 



most probably originate from the disordered component 
of the magnetic moments. The dual nature of the ground 
state is confirmed by the magnetic contribution c mag (T) 
to specific- heat [l3|, which is proportional to T 3 below 
TV and exhibits an additional broad hump around 3 K 
[Fig. Etc)]. The former implies three-dimensional antifcr- 
romagnctic gapless excitations [HI. [25l427j . while the latter 
reveals additional thermally activated exp(— e a /T) term 
due to the gapped magnetic excitations [261 ]. 

Finally, we point out similar spin dynamics in the mag- 
netic ground state of volborthite [22| , where spin-density- 
wave-like modulation 22j and IC spin correlations [28[ 
were found. Moreover, recent calculations [29[ showed 
that volborthite, like FeTe205Br has to be treated 
as a frustrated coupled-spin-chain system. This suggests 
that such systems are keen to form IC AMOD magnetic 
ground state, which offers a sound phcnomcnological ex- 
planation of the coexistence of PSD and LRO. 

In conclusion, we have found that in FeTe2 05Br the 
elongated-elliptical IC AMOD magnetic structure per- 
sists down to T/Tn ~ 1/200, as a result of frustrated 
chain topology. In this state fluctuations of the disor- 
dered spin component at each magnetic site are intrin- 
sic - they manifest as PSD at T— 5-0. Similar observa- 
tions in volborthite suggest that IC AMOD magnetic 
order is a natural habitat for PSD and can serve as a 
model structure inherently encompassing the intriguing 
coexistence of PSD and LRO. This conjecture could be 
tested on linarite, for which frustrated chains were re- 
cently reported to induce a LRO AMOD structure [3pj |. 
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Further in-depth theoretical investigations are required 
to account for PSD in IC AMOD structures on the mi- 
croscopic level. 
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NEUTRON EXPERIMENTS AND MAGNETIC 
STRUCTURE REFINEMENT 



Spherical neutron polarimetry experiment was per- 
formed for three different crystal orientations. In ad- 
dition to the hkO orientation (01), where the scattering 
plane was defined by the (100) and (0 1 0) reciprocal vec- 
tors [l[, the crystal was rotated to the scattering plane 
defined by the (010) and (10 2) vectors, (orientation 2 - 
02). Two additional orientations (03a) and (03b) cor- 
responded to the scattering plane with the (0 01) and ei- 
ther q a = (0.5 0.463 0) or q b = (0.5 0.537 0) vectors. This 
was achieved by tilting the crystal away from (1 1 0) by 
±1.7°. The 54 accumulated polarization matrices were 
obtained for 45 different magnetic reflections, i.e., some 
of them were measured in two different crystal orienta- 
tions or with reversed incoming neutron polarization. 

To estimate the experimental uncertainty and to elim- 
inate its influence on the modeling, we proceeded as fol- 
lows. First, we took into account leakage of the oppo- 
site polarization in the incoming and scattered beams, by 
considering imperfect, 97.5(5) %, efficiency of the instru- 
ments polarizers. Second, since synchrotron data [3, 0] 
showed that the change of the crystal structure at Tjy 
is very subtle and is not yet fully determined, we had 
to neglect the existence of the nuclear-magnetic interfer- 
ence term and to assign purely magnetic origin to all IC 
reflections. To avoid unphysical results, we therefore in- 
creased the estimated standard deviations of the matrix 
elements so that they within the error bar obey the gen- 
eral symmetry of the polarization matrix for the purely 
magnetic case, e.g., P yy = -P zz , P yx = P zx , ... 041]. 

Applying the above approximations we performed sev- 
eral simulated annealing runs, refining all datasets (in- 
cluding 87 integrated magnetic-peak intensities from the 
conventional neutron diffraction experiment) to the same 
magnetic structure model, but at the same time allowing 
different domain populations for each experiment, i.e., 
separately for 01, 02, 03 and for integrated intensi- 
ties Q. The minimized total cost function is defined 




as costtot=Y l jXj/(N johs -N pal ). Here j is the num- 
ber of the datasets and N pai is the number of the fit- 
ting parameters. For each dataset with Nj ^ s obser- 



FIG. 5: Pix (i = x, y, z) elements for [h k ± 5 1) reflections for 
+ and - incoming polarizations, signifying the effect of the 
nuclear magnetic interference. 



vations we define x] = Lil\° bs obs - X icaXc ) 2 / a 2 Xi ^, 
where X denotes polarization matrix elements P and/or 
integrated intensities /, and (Jx ioba is the estimated stan- 
dard deviation of the observation. For the best mag- 
netic structure model (Table 1 in the main text) the 
resulting total cost is cost to t = 33.7 and partial costs 
costp=Xp/(A r Pobs _ -^par) =9.6 for polarization matrices 
and costi=X//(-/V7obs-iVp ar ) — 24.1 for integrated inten- 
sity. This is significantly better than cost to t = 36.1, ob- 
tained for the collinear AMOD model with costp = 14.7 
and costi = 21.4 and that of the cycloidal (full moment) 
model costtot = 36.2 with costp = 14.9 and costi = 21.3, 
especially when considering that polarization matrices 
hold more detailed directional information. A small devi- 
ation between the calculated and the observed polariza- 
tion matrix elements is largely due to a weak nuclear- 
magnetic-interference term, which is most pronounced 
in Pi x , i = y,z matrix elements, where the first and 
the second indexes indicate the polarization of the in- 
coming and outgoing neutron polarizations, respectively. 
In particular, we where able to model the change of 
the Pix elements on reversal of the incident polariza- 
tion (Fig. [5]) only when considering additional nuclear- 
magnetic-interference term - the magnetic chiral con- 
tribution alone was not sufficient (3j-|6j. The apparent 
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nuclear-magnetic interference suggests the presence of 
higher order (> 2) terms in the magnetoelectric coupling 
that may explain a small induced electric polarization [l[ ■ 



fj.SH DATA ANALYSIS 

The use of powder samples is essential to ensure that 
exactly 1/3 of muons have initial polarization parallel to a 
local magnetic field, which is most important for further 
analysis. To ensure the quality of the powder samples 
several well-characterized single crystals were grinded. 
We took special care of the background signal, whose 
presence is inevitable in the MUSR ISIS (MUSR instru- 
ment at the ISIS facility, Rutherford Appleton Labora- 
tory, United Kingdom) setups due to muons stopping in 
the sample holder, cryostat tail, etc. As we combined 
data from two ISIS setups [a dilution refrigerator (DR) 
and a regular cryostat (RC)] as well as previous data 
from PSI, all with different backgrounds, these need to 
be properly subtracted from the data. The value of the 
background for each setup is set by two constraints: (i) 
from a non-frozen fraction at low-temperature [amplitude 
of the remaining oscillating part in a weak transverse- 
field (TF) experiment below Try - contribution of muons 
that are not exposed to the static electronic magnetic 
field from the LRO and precess around the applied TF 
field] , and (ii) by matching the ZF depolarization curves 
measured at the same temperature. We got about 15 % 
(depending on the setup) of background, which is typical 
for the used ISIS instrument (with given slit opening and 
sample size) and matches well with the known relaxation 
of a silver sample holder. Additionally, the subtracted 
ISIS and PSI data, all taken at T = 2 K, match well with 
each other [Fig.[6ta)] . 

Nuclear-magnetic field distributions at electrostatic- 
potential minima (Table lTl| were calculated by randomly 
orienting all nuclear magnetic moments (with appropri- 
ate abundance) within a fictitious Lorentz sphere with 
radius large enough to ensure convergence. 

Similar approach was used for calculating the elec- 
tronic magnetic field distributions at Pi and Pi sites in 
the LRO magnetic state. These distributions were ex- 
actly determined by magnetic structure model and ac- 
quire finite widths due to IC modulation of the magnetic 
structure along the 6-axis. The (contact) hypcrfine inter- 
action between muons and electrons is neglected. There 
are several justifications for this. First, both known muon 
stopping sites P\ and P2 lie well outside the exchange 
paths (outside be planes defining the plane of coupled 
chains) . The distance of both to the nearest iron is more 
than 2.5 A, whereas the nearest oxygen is 1.35A away 
(TablelLT| . Therefore, we do not expect significant elec- 
tron density on these sites. The latter distance is in 
reasonable agreement with the usual distance ~1.0 A, 
expected when muon bond to oxygen Q. Second, the 
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FIG. 6: (a) Comparison of the ZF fiSK data for 2 K obtained 
at different instruments and experimental setups after back- 
ground subtraction, (b) Comparison of muon depolarization 
in ZF and in longitudinal applied field (LF) of 1 mT measured 
at ISIS. 



calculated dipolar fields match perfectly with observed 
oscillations in muon polarization (with no adjustable pa- 
rameters). 

In search for the possible third muon stopping site, we 
calculated nuclear-magnetic field distributions through- 
out the unit cell. This allowed us to identify posi- 
tions with fields corresponding to experimentally deter- 
mined value. Since these positions are far from the local 
electrostatic-potential minima, we speculate that the cor- 
responding muons may deform their local environment 
and change the electrostatic potential. Given the high 
field strength at this sites, their positions have to be away 
from stopping sites Pi and P 2 in the vicinity of oxygens. 
One of the possibilities is that they are close to Br atoms, 
which possess the highest nuclear magnetic moments. 

Finally, field decoupling at lowest temperatures im- 
plies that relaxation of the "1/3 tail" is not due to 
nuclear field relaxation on the unknown stopping site 
^3/7^ = 0.304(5) mT] and thus show that it is indeed 
dynamical [Fig.[6][b)]. 



DENSITY-FUNCTIONAL-THEORY 
CALCULATIONS 

The electrostatic potential was calculated ab initio 
within the framework of the density-functional theory 
and the local-density approximation (LDA) Q of the 
exchange-correlation potential. We applied the Wien97 
code [10(, which adopts the full-potential linearized- 
augmented-plane- waves (FLAPW) method [llj. The 
calculations were performed for the structure with the 
experimental atomic positions and the lattice parame- 
ters, whereas the muffin-tin radia were 2.1 a. u. for the 
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TABLE II: The most pronounced local electrostatic-potential 
minima Pi, with corresponding positions, distances to the 
closest oxygen (do-Pi) and iron (dpe-pj atom, and local 
nuclear-field distribution widths (A/7 M ). 



i 


X 


V 


z 


rfo-p,(A) 


d Fc - Pz (k) 


A/ 7M (mT) 


1 


0.59 


0.13 


0.43 


1.36 


4.50 


0.071 


2 


0.25 


0.15 


0.90 


1.34 


2.55 


0.065 


3 


1.00 


0.13 


0.54 


1.50 


2.88 


0.017 


4 


0.91 


0.15 


0.66 


1.51 


2.88 


0.018 


5 


0.91 


0.09 


0.20 


1.42 


2.28 


0.039 



Fe atoms, 1.98 a. u. for the Te atoms, 1.5 a. u. for the O 
atoms, and 2.68 a. u. for the Br atoms. The plane-wave- 
cxpansion cut-off energy was set to 16 Ry, and the sum- 
mation of 333 k- vectors from the full Brillouin-zone (BZ) 
was carried out by means of the Gaussian method [12J 
with the smearing parameter of 0.02 Ry. 

The potential was calculated on a discrete mesh con- 
taining 100 x 100 x 100 points within the unit cell. 
From the calculated data we determined the most pro- 
nounced local electrostatic-potential minima given in Ta- 
blc|TTJ These were all found in vicinity of oxygen sites, 
as expected for oxides We note that potential wells 
around the listed minima significantly differ among each 
other, i.e., some being very sharp while others very broad 
[Fig. 3(b) in the main text]. Latter thus allow muons 
more freedom when choosing their actual stopping site. 



SPECIFIC HEAT - DETERMINATION OF THE 
LATTICE CONTRIBUTION 
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FIG. 7: Specific heat measurements (solid squares), model of 
phonon/lattice contribution (solid line), and derived magnetic 
contribution (open circles). Additional dotted and dashed 
lines indicate Debye and Einstein contributions, respectively. 



sponding four Einstein contributions |14| . As a result we 
obtained a very good description of the high-temperature 
data (Fig. [Il) with the Debye temperature Td = 280K. 
We note that low-temperature (T < Tjv) part of the spe- 
cific heat results can be fitted also with the Debye model 
alone, which would suggest that c mag (T < T/v) is ab- 
sent. However, such model yields completely unphysical 
behavior above (T > Tjv), making this model unrealistic. 



To deduce magnetic contribution to the specific heat 
c ma g(T), one needs first to estimate phonon contribu- 
tion due to crystal lattice vibrations. Since there exists 
no isostructural nonmagnetic compound for FeTc205Br, 
we need to model the high-temperature data, i.e., above 
50 K, where magnetic correlations should already be neg- 
ligible 13j. The most common approach is based on the 
Debye model [3], where vibrational modes are approx- 
imated by elastic vibrations of an isotropic continuous 
body, with a linear dispersion relation. However, this 
model alone cannot explain our experiments since to- 
tal specific heat (c p ) at high temperatures shows much 
stronger temperature dependence. In order to explain 
high-temperature phonon contribution we refer here to 
far-infrared investigation of the isostructural magnetic 
compound FeTe2 05Cl [l5|, where a number of vibra- 
tional modes were identified. The fact that most of these 
modes are sharp suggests that they correspond to well- 
defined nondispersive optical modes. Therefore, in addi- 
tion to the Debye contribution accounting for (linearly) 
dispersive acoustic vibrational modes, we took into ac- 
count the strongest four optical modes, i.e., 320, 400, 
470, and 650 cm" 1 , and approximate them with corre- 
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